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Preface 


In  this  thesis  the  methods  of  tensor  analysis  have 
been  used  to  obtain  expressions  for  the  basic  equations 
of  fluid  mechanics  in  terms  of  several  orthogonal  curvi¬ 
linear  coordinate  systems.  V/hile  these  methods  can  be 
used  with  any  valid  coor'^inate  transformations,  the 
reader  should  be  cautioned  that  the  specific  results 
obtained  here  are  applicable  only  to  the  particular 
transformations  which  are  listed  in  Appendix  A.  There 
are  other  methods  of  defining  the  various  coordinate  sys¬ 
tems,  but,  to  make  use  of  these  differing  definitions, 
the  reader  would  have  to  start  with  the  basic  equations 
listed  in  this  report  and  derive  his  own  final  results. 

It  should  also  be  emphasized  that  this  report  does 
not  attempt  to  explain  the  mechanics  of  tensor  analysis. 

If  the  reader  is  unfamiliar  with  this  branch  of  mathematics 
and  wishes  to  gain  the  background  to  enable  him  to  fill  in 
the  steps  which  have  been  omitted  in  the  development  of  the 
relationships  used  in  this  report,  the  book  by  Sokolniicoff 
(Ref  4)  is  an  excellent  text. 
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density 

tensor  component  of  the  velocity 
physical  component  of  the  velocity 
metric  tensor,  ^  ^  ^  ^  j 
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time 
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Where  tensor  quantities  are  defined,  the  symbols  are  the 
same  for  covariant  or  contravariant  tensors. 


vil 


GAH:/Mii:/62-4 


Abstract 

The  applied  science  of  fluid  mechanics  makes  use  of 
three  basic  equations  to  analyze  and  predict  the  state  of 
a  fluid  in  motion.  These  equations  are  the  equations  of 
motion,  energy,  and  continuity.  In  most  flow  problems, 
these  equations  cannot  be  solved,  or  are  very  difficult 
to  solve,  unless  they  are  expressed  in  terms  of  a  coordi¬ 
nate  system  which  conforms  to  the  surface  of  the  duct  or 
body  which  shapes  the  flow. 

This  report  utilizes  the  methods  of  tensor  analysis 
to  transform  the  basic  equations  from  their  Cartesian  forms 
to  expressions  in  ten  orthogonal  curvilinear  coordinate 
systems.  The  derivation  process  is  outlined,  and  the  final 
results  are  tabulated  for  each  of  the  coordinate  systems. 
Although  this  report  assumes  a  Newtonian  fluid  model,  the 
viscous  stress  components  are  listed  separately  so  that, 
given  the  proper  expressions  for  the  viscous  stress  com¬ 
ponents,  the  results  may  also  be  applied  to  a  non-Newton¬ 
ian  fluid. 
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THE  EQUATIONS  OF  FLUID  MECHANICS 
EXPRESSED  IN  CURVILINEAR  COORDINATES 

I.  Introduction 

Fluid  mechanics  is  defined  as  the  applied  Science 
which  deals  with  the  principles  of  both  gaseous  and  liquid 
flow.  Practically,  we  are  concerned  with  fluid  flow  over 
solid  bodies  or  through  various  types  of  ducts  or  channels. 
In  order  to  predict  and  describe  such  flow,  we  make  use  of 
three  basic  equations:  the  continuity  equation,  the  energy 
equation,  and  the  equation  of  motion.  These  equations  are 
dependent  upon  the  use  of  some  three-dimensional  coordinate 
system  to  describe  the  properties  and  the  movement  of  a 
fluid. 

For  simple  flow  patterns,  such  as  flow  over  a  flat 
plate,  rectangular  Cartesian  coordinates  are  adequate  for 
complete  description  of  the  flow  characteristics.  However, 
for  flow  around  more  complex  body  shapes,  we  find  that  the 
equations  cannot  be  solved,  or  are  extremely  difficult  to 
solve,  unless  they  are  expressed  in  terms  of  a  coordinate 
system  which  is  compatible  with  the  geometry  of  the  body. 

In  particular,  we  must  have  a  coordinate  system  which  has 
a  coordinate  surface  closely  approximating  the  shape  of 
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the  body.  This  is  especially  important  in  enabling  one 
to  express  the  boundary  conditions  in  a  simple  form.  Some 
examples  of  flow  situations  which  require  more  sophisti¬ 
cated  coordinate  systems  are:  flow  over  submarine  hulls, 
whose  shapes  are  elongated  spheroids;  flow  through  ellip¬ 
tical  pipes;  and  flow  through  converging-diverging  nozzles, 
in  which  the  nozzle  walls  approximate  hyperboloids  of  one 
sheet . 

The  problem  is  one  of  expressing  the  familiar  Cartes¬ 
ian  forma  of  the  basic  equations  In  terms  of  the  various 
curvilinear  coordinate  systems.  To  do  this,  the  methods 
of  tensor  analysis  will  be  used.  The  general  approach  has 
been  outlined  for  the  continuity  and  motion  equations  by 
McConnell  (Ref  1:  271-313) >  among  others,  and  results  have 

been  obtained  for  cylindrical  and  spherical  coordinates. 

The  methods  used  in  this  report  differ  only  slightly  from 
those  used  by  McConnell. 

The  fluid  model  is  assumed  to  be  viscous,  heat  conduc¬ 
ting,  and  isotropic.  Chemical,  electromagnetic,  radiation, 
and  diffusion  effects  are  Ignored.  The  fluid  is  also 
assumed  to  be  Newtonian.  However,  the  expressions  for  the 
viscous  stress  tensor  are  listed  separately  so  that,  given 
the  proper  expressions  for  the  viscous  stress  tensor,  the 
results  may  also  be  applied  to  non-Newtonian  fluids. 


2 


GAE/Mi]:/62-4 


II.  Basic  Equations 

The  basic  equations  are  presented  here  In  fairly  con¬ 
ventional  form,  and  they  are  then  adapted  to  a  form  more 
suitable  for  final  tabulation.  If  the  reader  Is  Interested 
In  the  derivation  of  the  equations,  he  Is  referred  to 
Appendix  B. 

Continuity  Equation 

The  continuity  equation  can  be  expressed  in  tensor 
notation  as 


where  the  i'th  component  of  the  velocity  vector,  and 
^  is  the  density.  The  second  term  represents  the  diverg¬ 
ence  of  a  vector,  and  It  can  be  expressed  more  conveniently 
by  using 

fV.)  ,  (2) 

where  V-  represents  the  physical  component  of  the  velocity 
vector.  Substituting  equation  (2)  into  equation  (l),  we 
obtain  the  continuity  equation  in  its  final  'form: 

3  P  i  £  /  /F  \  /  \ 

??  [f^,  fVi)  =  0  .  (3) 
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Energy  Equation 

In  tensor  notation,  one  form  of  the  energy  equation  is 

j>[ft(k'  t  Z  y;i7‘|j]  =  If  * 


where  h  is  enthalpy,  if*' is  the  ij'th  component  of  the 
viscous  stress  tensor,  p  is  pressure,  and  is  the  i'th 
component  of  the  heat  flux  vector.  By  substituting  physi¬ 
cal  components  and  making  use  of  the  general  form  of  equa¬ 
tion  (2),  we  obtain  the  energy  equation  in  its  final  form: 

+  Jx‘6l  +  i  -  ?t 


-L  9  .  (  /T  A  1/  p  1  d  I  ff  \ 

Ktr  V ,  (5) 


/T 

whore  ^ , 

Fi. 

whore  It  \  ,  and  are  the  physical  components  of  the 

viscous  stress  tensor,  the  body  force,  and  the  heat  flux 
vector,  respectively.  The  ij'th  component  of  the  viscous 
stress  tensor  can  be  written,  after  Sokolnikoff  (Ref  4: 
321-324),  as 
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5  1^) 


'T' 


■f 


> 


(6) 


where  Is  the  coefficient  of  bulk  viscosity  (usually 
taken  uo  be  zero),  and  ^  is  the  coefficient  of  shear  vis¬ 
cosity.  After  carrying  out  the  indicated  operations  and 
substituting  the  physical  components  for  the  velocity 
vectors,  we  obtain 


=  (■j- It) 


2  LJ.‘.A 
/ 


(  c.i 
I  q  ^x^ 


1  (  tM  j/  Oyw  t4\  _  ? 

3  M  a  'ipSTJ  • 


(7) 


V/e  can  obtain  the  physical  components  of  tr^^by  multiplying 
equation  (7)  by  •  Although  the  multiplication  will 

not  be  shown  here,  the  physical  components,  rather  than  the 
tensor  components,  will  be  tabulated  in  the  next  section. 

Equation  of  Motion 

The  tensor  form  of  the  equation  of  motion  Is 
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The  I'th  component  of  the  acceleration  vector  is 


1  dy;,i> 


(13) 


The  final  form  of  the  acceleration  vector  is  obtained  by 
expanding  equation  (13)  and  substituting  the  physical  com¬ 
ponents  of  the  velocity  vectors: 


The  term  on  the  right  side  of  equation  (12)  can  be  ex¬ 
panded,  using  equation  (9),  as  follows: 


where 

1,'.!  h" .  W  . 


(16) 


Expanding  equation  (l6)  to  eliminate  the  Christoff el  sym- 
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bols,  we  obtain,  after  simplifying 


i  3X‘ . 


Substituting  equation  (l?)  into  equation  (15),  '•;c  obtain 

2.  »  J  .  ] 


Substitution  of  the  physical  components  of  the  '/Iscous 
stress  tensor  yields 

V  -  --L  f.n 


Ali'< 

5X'"  . 
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III.  The  Completed  Transformations 

For  convenience,  the  equations  that  are  to  be  expres¬ 
sed  In  the  various  coordinate  systems  have  been  restated 
on  a  fold-out  sheet  In  Appendix  C. 

Cartesian  Coordinates 

x'=  x  x'=  4  x*=  Z  (so) 

<3.,-  =  <3  =  (21) 

Continuity  Equation. 

It?  *  li(jV.VaVpV'3)^l(fV.)  -  0  (S3) 

Energy  Equation.  The  left  side  of  equation  (5) 
becomes 

J  I  li  + 1  (v/+  +  Vx  Fx  [  +  I  ( ^4^  + 

i(^x  ^  V'j’’  +  v'^'’)]  \ 


+ 

(23) 
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The  remaining  terms  of  interest  arc 


^  VjV'-w  V.fr-'O  + 

*  V,P‘)  ,  (24) 

^V.hi-  ^  (v,  F,  +  +  Vjf,)  ,  (25) 

and 

^'l  =3X  (pt  |j(^,V|j(p  .  (26) 

The  components  of  the  viscous  stress  tensor  are 


tr''-- 

-- 

/3V^  ^ 

'^\9X  ■^^ 

SVx  \ 
ajj)  , 

(27) 

tr"'  -- 

-- 

\  ?x 

Pi  ;  , 

(28) 

= 

p?  j  , 

(29) 
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=  (>2-  3^)1  9X  ^  + 


ii'*  ]  .  7  ^  ii'- 
5i  J  ■*  ^'r 


/  i  s  [  ay.  sv; 

('?-  jT)  [  5X  +  + 

^•1  X  ,  5Vs 

Sii  + 


^Si  -  /  ^  ^\\  2y*  9V^ 

^  -(>j>-3t)[^x  ■'‘a^ 

Si  J  *  .  (33 

Eq^uatlon  of  jtotlon.  The  components  of  the  viscous 

stress  tensor  are  tabulated  above.  The  remaining  terms  of 
interest  are 

K  -  ^  V.  ^  1^^  +  y^  g  ^ 

li  '  3/  K<  Ja  Vi  +  v;  1^*  ^  (33, 
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4  * 

7^.  \i..  -  -  ay  ^  ^ 


Cyllndrica 1  Coordinates 


X 


I 


V- 


-  £ 


S"  '  ^33 

q--  =  ? 


1 

7- 

X 


Continuity  Equation. 


GAi/NV62-4 

ft  *  i  lk(ff  V,)  •»  r  (f  Vi)  =  0  (42) 

Energy  Equation.  The  left  side  of  equation  (5) 
becomes 

f  {st  l.Vi -t  V,‘)]  + 

V,  fi  [li  +  i  (v^'t v/t  (i')] I  .  (43) 

Also, 

fx'  V,'  ?r‘0  =  f  It  [  >■  ( Vi  IV"-*  ir’% 

Vj  IT  *‘')]  +  7  f'e  Ufr"*  V,ff*V  + 

,  (44) 

7,  =  tVsTs  +  Vjf,  )  ,  (45) 

and 


13 


GAE/MS/62-4 


The  comnonents  of  the  stress  tensor  are 


fr 


*'*  -  tr  ^2•  J  , 


f;:®*,  -)?.*»  .  v[t  l^e  *  H]  , 


tr 


yir 


J 


^  1  £y« 


0  3Vr 
+  2.T 


'  C'^-  |r)  [  y-  a)'(i"V^)  +  H 

iYil  o  i  \/  1 

^  Irl'y  Te  ^Tt\  , 


(47) 


(48) 


(49) 


(50) 


(51) 


and 
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^  MoUon.  In  addition  to  the  components  of 
the  stress  tensor  tabulated  above,  the  terms  of  interest 
are 


A 

It 

vi^.Qj 

*  K 

i  U  ,/ 

1-  *«  5e  + 

V/ 

“  T  , 

(53) 

^e 

■■  H 

1/ 

'  h  3^  , 

(54) 

A 

3t 

+  \/ 

is  ^  , 
f  5e  V 

f 

z  j 

(55) 

£ 

-  -it . 

)  +  Tc^e(^> 

yo\ 

) 

-  h  ir®® 

i 

(56) 

QiK 

ik 

■■ 

i 

"  "  V" 

1  +■ 

k 

> 

(57) 

and 
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Jherical 


Coordinates 


'5"  “  I 

^33-- 

^  =  1^''  Sin'^e 


(59) 

(60) 

(61) 

(62) 
(63) 
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Enort^y  Equation.  Tho  left  side  of  equation  (5) 


becomes 


j  {at  [  h  1 1  ( v/*  y'j’’)]  +  K  li-  [  In  ♦  X  ( v',.*’  <■ 

V*  *  p*  re  ^  I  (v.%  V*  i/,-)]  f 


y? 

Y'  S/n  © 


Also, 


i  IJIL  t/  /^  ••\  9  r  X  /  A 

<5  axHq;];i]7  Vetr^j)  ^ 

FsTiTe  ^ 

"■  fhrs  ,(66) 

■j*  ^  '  f  ^  ^1-  ^  ■*■  fvf  J  ,  (67) 


JL  ^  /^\  A^/x-An 


|r  3IM  6  ae  (sme  ^  0  ^  ^ 
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The  components  of  the  viscous  stress  tensor  are 


=•!>«'«  =  T  [  ^6  (  57^  )  * 


-i-  ^v»-| 

V  5fH  6  J  J  (71) 

(^- It)  It*' Jh(l"*'VJ.)  ^  ’TfiKliel^Q^inB)-^ 

~hn  e  ijf^]  *  ^  (78) 

(i^-  It)  [t'-  li-Cif’'/)  ■*ir^9lek^"’®)  + 


_l _ 

V”  J/n  6 


_A_  .  _  r  i- 


Y  Sin  6 


5i<pi  r  i  /■5\4, 

3?  J  *  FltTfl  (3/  + 

v;.  ■>■  14  c.o<i  fl)]  ,  I 
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Equation  of  Motion.  The  terms  of  Interest,  In  addi¬ 
tion  to  the  components  of  the  viscous  stress  tensor  already 
listed,  are 


f 


„  ay.  ^  0-. 


i- 


%  0.  _ 

r  Sin  6  y*  j 

a'i  *y;-a  /  ay*  ^ 

•e  dt  r  dr  irvj*--p:  gg  + 

V  d9(^in  9  I/,) 


(75) 


(76) 


v,  3y» 

Sm  6  9^  j 


(77) 


.JX 


Tii.K  =  -ir  +  r' dr  ( r'"{r''')  * 


r  sme  90  (^i>^  0  ^r"®)  + 


(73) 
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1^C 


Ui,K  =  -  ^  30  Sr  irk  J'- 

I-  Sin  s  3«  ®  '•'  fsirTedilk^'^)  ~ 


and 


Cot  0  +  Z  )?■ 


6k 


(79) 


-  3  *  1  Op  ^  /  . A  _  \ 

^3i.K  ""  ~  i'  driir'k^  )  ^ 


vleiP^) 


1  d 


Ir  Sin  e 


Z  CCJ  0 


(SO) 


Parabolic  Coordinates 


x^= 

X^=  (p 

(Bl) 

>  •»  M 

(B2) 

A  -f  A 

(83) 

^33  =  Aa( 

m 
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(At 

— 


Continuity  Equation. 

ST  *'■  ^Tm  aA  f  Va]  + 

ff  =  0 

Snorgy  Equation.  The  left  side  of  equation  (5) 
becomes 

flttlu.  i(v,\  C 1/;)] 


Also, 
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.  (93) 

=  (if  -  |T)[,r^  Vi)  + 

'rf'f  )  ^ 

'  .1  — J:.3^__ .  1/  7 

V  ■  /r(Af^)''’-  ViJ  ,  (94) 

r  iL  9 

=  (^^-  |t)1  dl  (fJJXTJ^  \(*)  + 

A^  +  jfe  l-H  +  iV 


%Jl 

(.a+m) 


i-'-V'  -  V^]  , 


fr'’’''’  =  (if  -  |r)lAv.  Ia  (/aTaT^  V^)  + 

_L  i-  /  , - i;  \  -i_  2y<>l  (-?:-  SY,  ^ 

A*>*  9/^  (i/m{Ai-/<)  \(«|  -I-  /Xm  ixp  J  +'l-LiCi7r  3^ 

VA(A+/()  Vx  ■*■  //^(A+/x)  .  (96) 
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uA-;, 

latl' .1  __f  i-:  11  ! , .  Km-  :  jmi  ont  iiii)  of  ch  •  vi.icous 
3lrei)3  tensor  are  listed  above,  'ine  remaining  terms  of 
Interest  are 


A 


If' 


,/  2Y*  wv  jx  ■■>  1 

^,7^  a<f  /A?^)  , 


(97) 


■-  at  xrjK  aA(tCwr IY^Ja^m 

1.  '•  /^  >  i 

fAM  «a<p  j  (93 


I  =  - 


f  (/XV^ ) 


v^kiorK)] 


Vo  ^v, 
+•  /Tm  3<p 


5  A  I  A  '  p  ^  ^ 

=  -2jFi7iaA  ■*-  /XfM^aAlAlr  7 


aIaoToVa.  ^M«_ 

OTaF^)  1^”’  •'•  ca"^)*'-  fr"*'  ,  (100) 


i  <? 


2i| 


GAE/MS/62-4 


gJK  _ _ 


Z  9 


(7^)’^-  ir^  -  - 

+  a^f'’  ,  ( 


(101) 


^  ^  _  —•^  -  — P  /X*  ^  /  A  ,v 

^3»  Ni,K.  “  AT''  d(p  ^  5A  (A^  fr^  )  -^ 


A-f 

■*■  .  (102) 


P^'olate  Spheroidal  Coordinates 


X  =  A 


X  =  f 


X  ^-h' 

In  =  ^  A  -  I 


a  =  cl’- 


(103) 


(104) 


(105) 
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'jAr,/HrV6i'-'i 


^33  “  6L  ^  A  '  t  )(  1  '  J 

(106) 

(107) 

Equation. 

al"  ^  l/oS7)(F^  pVj  + 

|a  l/wvxi^  f  V^]  + 


-  r-  ^  £ 

aip 


[pV'f]  = 


(108) 


aquation .  The  left  side  of  equation  (5) 


becomes 


platih  *  I  (y;%  v/a  +  cL  Z^'-  ax  [a+  i  (y;’’+ 


"")]  +  a  |;  [  A  *  i  (//*  v/f  K/). 


I  ^v;}h^U 


CkfCA 


-h  +  Z  (V'^^+  O-  ^  • 


(^09) 


Also, 
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Th6  componGnts  of  thG  viscous  strGss  tensor  are 


3  /  i 
5A 


(113) 
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and 


2*  jEI,  w  15  ^'■■1 


(U8) 


Equation  of  Motion.  In  addition  to  the  components  of 
the  viscous  stress  tensor  listed  above,  the  terms  of  Inter¬ 
est  are 


29 


Vl  -  l*i  , 


,  y,  io:  ±  _ _  . 

dt  CL  +• 

^  V</-  9V^ 

CL  ih^'-i  +  oTfAMHi^  d(f  ^ 

iiYji"  -/nr^ .  mc 


P'  syip  9  / _ _  ^ 

\  ■-  ^t  ^A  Vcp )  ^' 

3^(/T^  \4)  -I* 


(120) 


a  v"(  A'  -jyi-  j 

-ll\  i  fA^  if _ 1  3  f  A  V .  1 

'ij^K  "  '^^AV^aA  l(A^-i)tr""J  + 


(121) 


ff*'')  +  If  (fr 


Xm''-i) 


>/a^ 


(-22) 
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___ 


M. 


lA/KT  A 


AtA 


(123) 


,  _ 1 _  3p 

Ui.K  ^/(A*'-i')a*V^  «3(:|)  •»■ 
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u.  ■ .  I'.i'./o 

Spheroidal  Coordinates  (Oblate  Spheroids) 


K'--  A  ^ 

(125) 

1m  =  A*--! 

(126) 

Su  =  l-M'' 

(127) 

^  y 

(128) 

C\'-i3a-Ai‘3 

(129) 

Continuity  Equation. 

3  p  fy- 1'  9 1  ^ _  \ 

at  ^  j  VJ  + 

sxjx  aif  ( j*  /(p)  ■  0 

(130) 

Energy  Equation .  The  left  side  of  equation  (5) 
becomes 
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Also, 

v;^wj^  .  (135) 

--  f(^A  ^y9^„)>  (133) 


and 
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The  components  of  the  stress  tensor  are 


3 

CL  3a- 


'(h)  + 


i  jy.  7 
^ A H  J  j 

[SiSflttf,) 


It"’  = 


-i-  ^  1 

CL\M  3(p  J  j 


(135) 


(136) 


(137) 
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-  |r )  f  /A'‘-/<‘  Vj, )  ^ 

A.  (  f- — ,  ,,  \  ^  9Vif  1 

9m  -m"'  v^)  -^oA/f  9(f  J  ■*■ 

r  ^  i  /Eir  \/  \  2^^  , 

t  ^  ^a(/a^-m‘  m  j  “  d  “ 

2_M  \/  1 

^  CA'-a'-)  '^  7<J  ;  (138) 

a  (a/P^  Va)  + 

|:t(M/A^  V>,)  +  aTJ 

(12/  fUST  „  \  jGEE,  ,/ 

L  a.  yHj  -I-  a  M  + 

^M  A  *  1  . .  1 

a  W  >  (139) 


and 
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^  |^r)[aAa^*)  3a(a/F^  1^) 

,,/r^  3-  /  , _ ../s  j  iy^i 


PTv  1 

r*^  v*  j  . 


(140) 


Equation  of  Motion.  The  Important  terms  In  the  motion 


equation  are 


f  ^  V,  nyrr  ^  ^ 


A"-V  9a<(/A-^'  yJ-^oJ^M  d<p 

Vx  jg^ ,  fiinr 

a,  olA  'Ix^'m'’  > 


4. 

■^- 


/O’!  (r- - -  U  \ 

CL  X^-M^  IfA^  K) 


^  IL^  iY^  ,  £Ym  ^ 
a.  d-H  aXu  a(p  ^ 

“OL  (A''-/,'-:)'^"  aM  /A^->q^  J 


(141) 


oT  (A''-vr"  aM 


(142) 
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r  11‘f  /A*-/  '  £  /  .  V 

-  9{  +  a.A  9a(  a  V,,,)  + 


(143) 


^  "hi.K  =  +  cuTJ^fxix/x^  ir^^)  i- 


,v  C  a/T^  .a 
ol /A*'-f  (A*  ,4‘)  (T  -  - 


A  ^  ^  ipii’ 


fl-A  /a’-'A'-  n  ^  ^^‘  a.  ^T*  , 

A.  i/HST^P  1/A^£/  a  \ 

•'S^  ^2i,K  "  ”  ^•VA-ah*'  5/<  +  axh^-M'-  ^A  (Afr^^j  + 


(144) 


— -^(A*'-1)  ^A  /H/IvH^  a  .  . 


1?^*^ 


uIEU, 

+  a  (A '•'A< 


(145: 
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Parabolic  Cylinder  Coordinates 


X 

il 

X' 

;(^--  2 

(1^7) 

X  +  M, 

= 

(148) 

A  •^/^ 

•(149) 

=  i 

(150) 

CA4./.3'' 

(15]) 

Continuity  Equation . 


^  gl/X*  ^  \ 

dt  X  +M  c?A  ( fXTJT  ^  )  "^ 

KTJ\  «9m  (/^TT  ^  1^) 
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Energy  Equation.  The  left  side  of  equation  (5) 


becomes 


y,')]  + 

ziy/t-y/i-Vf")] -f-  zy^/^iz  ^  [h  + 


(153) 


Also, 


Vm  '‘V  li  )t  '•‘)]  >7^l<  [/;n^  fi/.  f?  4-  vL 


-»-  Ktr^-h 


f  K  '-'.  =  -f  +  ^-  V'i.?,)  , 


(154) 


(155) 


A  £.  /  /T 

X/m  5  /, 


f)  = 


a/7  2 

/  •+/U  <?A 


(156) 
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The  components  of  the  stress  tensor  are 


A  A 


'[2/T9A{f 


i  /-JL. 

3A  ( /JtXT 


It  ’  ?T 


(157) 


(153) 


r,  HET^v, 
l^htM  9m  J  , 


(159) 


1T**=  |T)f^lA  Va)  + 


j-pMM 


Va  +■  , 

^  (>?-  3T)[a^ 


(160) 


IdR 

StM 


i^)+ M'J  +  ^h§JI^  1^) 


7“  fX  II  *A,A  1/ 


(161) 


4o 


tr""  -  (>e-  ir)[^  Ia(/a^7  Vj  + 


1^)  +  J 


Ji  . 


(162) 


equation  of  Motion.  The  important  terms  in  the 
motion  equation  are 

? .  iy-*  _ ,  „  or  ii,  9 , 

>  <?f  9a  9/(  4 


32  -Vm  oOi*'’-  , 
^  ^  V  iZ  S.h _ _  \ 

3t  ^  M  A+4(  3A  ^Ai-M  \^ )  + 


(163) 


T  /  /TIO  =r^ 


(164) 


G/Aiv'I'i2/n2-4 


/§.  sJ/T 'h-^) 


-S- 


Z/XT  A 


(AtM)^‘'  -h  aTTo*'--  h*  ^ 

7^”  *•  z/j^ fA(lr'“'')  + 


(166) 


7Thk((p:  Ir'^'')  ^  1> 


2  /A>?\ 


ifM'dtM')*'*- u^)^'-  1t"*  , 


(16?) 


/fc  Vk=  -li  ^  X^Ia(/h71t«')  + 


(168) 


Elliptic  Cylinder  Coordinates 


X  --  A 


X --  e 


aL'‘-(A^-Mp 
%--  X--1 


(169) 


(170) 
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.■•s^A.'Vf)2-4 


a‘(A^-.L4*) 

i-M"- 

(171) 

1 

(172) 

CA‘-i)(i-^*) 

(173) 

Continuity  Equation. 


i?  ^  /iFT  a 

at  CLCA'-^M  5>> 
7 1  -  ,u^  3 


(/F^  f  Va) 


+ 


Snercjy  Equation.  The  left  side  of  equation  (5) 
becomes 

[It  [ lo  ^ V^‘+  Y,')]  +  a  /^.  Ia  [  + 

i(VAV\^'ty,')]+  [  u  + 

Vi)]  +■  az  Jn  1 1  (v^%  ^ 


(m) 


(n5) 


Also, 
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The  components  of  the  stress  tensor  are 


(179) 

(180) 

(iSl) 
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dH  {(T^  l/^)  +!?']+  l/J  - 


-JtiiUjiMiL-,  w  iiiflsr,  ,/ 1 


(182) 


|r)[a.(A^‘)^A(/A^  \4) -f 


(133) 


and 


tr**’  (>^  -  |r)I K)  + 

£  I’n — , ,, \  pi4 1  3Vi 

a(AV'  a/<WV  si  \^  . 


(1S4) 


Equation  of  Motion.  In  addition  to  the  components  of 
the  stress  tensor  listed  above,  the  Important  terms  In  the 
equation  of  motion  are 
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ta  =  at  ■'‘a. ^  +  o:  "a 

w  ^  jEET 

14  3?  ~  a  a*x4"‘  1  (135) 


GAE/ME/62-4 

and 


(190) 


Ellipsoidal  Coordinates 


X'i  ^  ^  y^z  j)  (191) 

g„=  ^(d.^-AMb'-AUt'--A^  (192) 

9x1=  ^  (193) 

(a^-yXb'-y>(c'-;^)  (194) 


9=  ^ ^  (!a.'-A)(b'-A)(c^-A)(a'-MH b‘'x<)(x<.cO(a.'-y)a'-y)^c*-v)(  195) 


Continuity  Equation. 
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i:lquatlon.  The  left  side  of  equation  (5) 


becomes 


p  [at  [  li  t  2  V  Vy'')]  + 


(197) 


Also, 
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A  —  ,  /  /  5  I  /  a''  1 1  \  ^  /  fA'*  A)^b‘'  >)/c  9  r .-  /  AA 

^<5xM/j;j;7  K  rr  V  =  ^-mVa-v)  ItMa-y)  I  V^7r  4- 

,A  .  a^i\7  2^(gL-->4Kb-->L<)U-(;*.)  9  r. _ 

Vu  tr^  )]'*'  (M-A)(y->L<)  <3/< L Y(h->‘)(^-A) (Y^It  + 


v;  1r  *’^)] 


i/^?r''%  v;  ir  '''Jj  +  (v-aSo.-:;)  -  k  v ^ 


(19S) 


(199) 


J:.  9  /  /T  \  l/{aM)(b^-A)(c«--A^  j.  /„ 

^5  ax‘  1/^73^  f  ‘  i  -  7A-/A)(A^y)  9A  l/c7 


vn)a-y) 


a.*-M')(b*'-AT(M-c^)  _9  /  _ _ _  A  \ 

(M-A)(y-/<)  i ^/H  j 

Z  1.  / _ _ _  A  \ 

y-A)(y-^)  9y  \  n>>->)ii>-M)  j  . 


(y-A)(y-^)  9y  W(y'A)(v-/M) 


(200) 


The  components  of  the  stress  tensor  are 


/jY  ^  ^  ^ 


[  /7a-.A)<b--A)(c^Ar  9  ^ 

IM  >  -  A  aA  I/mTT  vIu  )  + 


j  ^  / _ 1  \/  \1 

/  7=77  (7^  YaIJ 


(201) 
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GA.-l/Mi/o?— 'l 


a''»A  r  (H'-AW-iHc'-M  w  ) 

Tr  =  ^  =  TUi/  <9^(7?^  lv/  + 


/ 


9  / _ 1 


a>» 


(202) 


r  fe  |/J  4 


Zv  av(?>^  j 


(203) 


/  2  J^^^<l^-AK6»-Ayc*-A)  9/ _ _  V 

It  -  ( >^-  I  ^)[  d^{fu-M)(A-)>i  y^j 


ZkQ.^-A)Lb^-M)(M-C'-)  9. 


-><)  ifiH->^'){y-M)  + 

'  ^  1_  /  /n‘'A)(bMKc‘-AV  ,  ,  \ 

^  aAi/  +  A 


Z/(^a.’--A)(b^-Ayc^-Ay  (2/l-A^-y) 
(y-A)*'*- 


V. 


Z  iio.'-x)(b^-\)  ^(&^->^)(c^-x')-^  Cb^^M(c^-x)} 
/(aJ'-y^)(6'-A)(c^-A)(M-A)(y-A3 


—  I  (^^■M)(li^-^)(c^-^)  , 

A 'A  v  (A->i)(y-/H)  ^ 


_i_  /(ci^-)>)(b^'y}(7^  ./ 
a^'A  f  (^.y)(M-yy  V 


gas/Mc:/62-4 


1  i  ((L"-A)(b^M)  ^  4.  ‘'-xv )  } 

\/ (CL^-m)C  ( y-zM^ 


_L  /(a^>yKb^-y)(7^  ,  1 

A->K  V  (A-y)(M-V)  V;;  J  ^  (205) 


and 
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tr  -  (>2"  ivU  ry-AUM-A)  5a (ir(M-A)(y*Ar  1^) 


2.  ^  /  N 

C  y->u.)(M-A')  <9m  ( l/  fyM-/<)(y-x)  ly  )  + 


2  /7a*-y)(6*-y)rc*-y)^  £.  f, _ _  ,,  \1 

VA-y)?M-y)  By  [u)f-A)(>f‘My  VyJl+ 

T  <3y  lif  a-y^M-yj  Vy/-+- 


_2_  (dL"-AVb*-AycM^  w 
X-y  1  /'u -A'i/y-AX  '> 


A-y  t  (M-AVy-A'i  '^A 

/  (Q.'‘-m)(6^-m)CC'->m)  sj 


At-y 


M  ■•“ 


lUo^-y)Cb^’y)Cc^'y)  (2y-A->u) 

(y-A)»/>-  (y-x)^-^^  ^y  ■** 

Z  [(a^-y)(b^>y)+  (a.^-y)(c".y)  c  6»-v)(c *-yj3  “i 

i((i^-yXb'^^y)(c^-yHx-y)(M-y)  .  (206) 


aquation  of  Motion .  The  Important  terms,  in  addition 
to  the  components  of  the  s'.rcss  tensor  listed  above,  arc 


52 


GAE/MS/62-4 


+*  =  (h-xw-T^  s>.  + 


7.  ((a}-yw-M)it^  ^  1 1 - - ./  \ 

A-  A  I  M - y  <5m  \im~^  1/^} 


/A- A  I  >M-y  uA^'A  M/'’’ 

M;  ISSUES' 


V  ^H-Ai(y-Ar  Y  (M-X)li>‘A)  3  (207) 

SVm  IVa  l(fl.^Ay b‘--A)fc‘-Aj'  _£  /  - - -  \y  \ 

T^-  di  A-m^I  y-A  aAl/Ai-A  VJ-»- 


<3  Va4 


^  ^  V  (A->m)(v'">H 


(f^  %)  - 


^-A  \  (A*>M)(y'>M)  "x-y  V  ,  (208) 


9t  x-y  V  a«-a  aA 


y-M 


f(a.^-y)(b^-v)(c'*--y) 
(A-y)(M-v)  ay 


X  ico.^-m'^-yKc^W  jC 

v-A  V  (A-y)(A<-y)  "’>'->M  V  (A-yyAt-A'i 


(209) 
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GAS/ME/c2-4 
5J1C 


A  A  _  _  -  j(&>Z0ib'^M)CC':Z)  i.p 

MLi.K  ^  t 


—  + 


vJm 


+  /p:T  U-aV'^  ^  ^- 


/7^  (y-Ai 


^  ?r 


(211) 


and 
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^  iji.K  (a-«cm-»  ai- 


Confocal  Parabolic  Coordinates 


A  =  A  =  y 

(213) 

(A-^)a-xj 

g,,  =  ^((l-A)Cb-A) 

(214) 

('m-aVm-v) 

(215) 
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(J33  =  -^(a-y)(b-vi  (216) 

g  '  4-^  (o.’^)(b~X)(M-6.)(b‘M)(Ci’y)(b- y)  (217) 


Continuity  Equation. 


■h 


^M->)(y-»  ( I'ryM-Aify-A)  f  + 


2jjM-a.){b-M)  i.  /  ^ ,  .  \ 

(A-M)(M-y  .3a  n(^->^H>u-y)  + 


z/(0i-y)(b-i^y  2  I _ _ _ 

(A-y)(>M-y)  <?y (va-yjfx-y) 


(218) 


Energy  Equation.  The  left  side  of  equation  (5) 
becomes 


f^at^^l  +  x(v';,^V^Vl//)]  +  IV^  [in  + 

itVi  +  K;  )]  +  ^V'y  i|(y-A)(y.^)  9y[  In  i- 


(219) 


Also, 
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BX 


<l/C>4-a)(b-A<^ 

TTv^TIjr^ 


U 


(\^fr '^+- 


./  ^  yu\l  Z/(T-yXb-yy  3^ 

+  V.tr  a-y)u-v)  3> 


^  ^  VyF,]  , 


(220) 


(221) 


1.  3  /  /F  ^  \  Zi/ca-Ayt)-A)  ^  ^ _  A  \ 

=  Cm-aMv-a)  ^a) 


Z  /(u-a.)(b-iO  ^  ^  \ 


U 


a-y)U-y)  \uA-v)(yM-)») 


(222) 


The  components  of  the  stress  tensor  are 


(M-(3^)Cb-AV  ^  (  —J:~ 
M-^  \  yA-aa’ 


(223) 
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j  (6.-y)(b“y: 

^  i  M-y 


lylThK)]  , 


(ri-  Ir'jf  7^^^  V^)  -h 


Z/rM-flL^^b-MV  3  /  - - -  ./  \ 

AnuVy-yK)  v)  + 


UTy)  Vj] 


Z^  U-y)(b-y)  i. 

CA-yVx-y) 


(224) 


(225) 


IaLIREMEI^  \/ 

T  aA  I T (A-M)(A-y)  ^A  +  U-mY'^  m  “ 


-:^(^A-d-b)  z  lleL-M)lt>-M) 


?===±===»«aipi— »  \/  C  {.CL-Mno-M)  X, 

^ra-A)(b'A)(A'>M)(A-VJ  V;^-  V  (/H-A)(/Ci-y)  Vm  “ 


K&-y)(b-yy  “I 

A-y  V  Cy-A)(y-/^)  Vyi  J  ^ 


(226) 
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h- 


Mm 


(a-aKA">») 


G(A-x)(A-V) 


f 


+ 


(M-y)(A-y) 


ayf/(A-y)(A4-v)’  V'y) 


•H 


[  3  /  f((l'M)(b-j^ 
tL^  a;:^W(M-AMM-y) 


4- 


Z  V(M-d.y  b-zMV  \/ 


VTd-MU  b->M)(yM-- 


Vu- 


^  /7Q.-y)(b»yT  ./ 

M-y  V  V'^  J  ) 


and 


(227) 


ga2/m::/62-^ 


ay^rzTTFgr 

^M-V)  (V-M*) 


BMiUx-M)lM-y)  ^m) 


+ 


Z^(L-v)(  b-y) 
Cy-x)()>-AA) 


(A-y)/>M-y) 


+ 


i-  / J<fl.-y)/b-yr 
3v  \i  ('y-x)(y-M) 


\U-M)ix-y) 


2.  fM^MjTh^  V/  Z/ra.-»rb-y)  (iy-A->4) 

Ym  ix-yy'^  yy- 

g  (g  y-d-fe^  \/  1 

Vffl.-y)a-v)fy-AUy'><)  V  J  .  (228) 


Equation  of  Motion.  Tho  terms  of  Interest  In  the 
equation  of  motion  are 


f  ?yA 


(a-A)(6-A)  ' 

5A 


Z  Vm.  I (M-dL)(b'M)  ^  (  I - ,  »/  \ 

X-M  v  V;^)  + 

f(i->>)(b-y)  d  (  , - -  w  \ 

A-  y  T/^-y  ay  WA-v  )  t 

ila-AMb-AY  ^  I  ((L’X)(b-x) ' 

A^-A  V(A*>M)(A-y)  y-A  /  (A-M)(A-y)  j 


(229) 
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y  Ra^yXb-y) 

p;^  )  (231) 
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^  l(a.-x)(6-A)'  £p  ^ 


ilf^w-AXb-A)  If"^)  * 

Ai=?Pfc(cr7f;>'). 

(lA  -M-vW(<X.-A)^b  -A^  ^ 


Q.-h  b-  ZA 

/f<i-AMb'AKA-M>a*vi  rr 
/^A-AiTb-A^  a2  >'V 


2.  licL-yXb-y) 
(A- )>)’'“  »~m7? 
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9^'*^ 


li,  K 


U&.-M)lb-M)  £p 
^  y  (m- A^M-Vi 


z.p 

"■ 


(A-M)(M-y)  'fr^)  ^ 


^  f(dr-x)ib-xy  d 

A-y  V 


A->m 


dxifA-/  1r^^)  + 


_L  Ra-yUb-yy  ^ 

y-A  Y-TTrjr-  (/r7 tr + 


6_  +  b  -  •Z.>cA 


7(67^!^) ^  yy  fcl-XUb-AVy^P^ 

Tr  -  fr  - 


.— ^.  ,,  ICa^yJ(b-yy^M^ 

(H-yy*-  I  Av- y  lr  ^  (233) 


and 
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hi,K  ■  "■^((>'-A)f>'->()  ay  + 

Wv-rt(y->4)  ay(/(a-y)(b-yy  ?T  **0  + 


/^a.-y)(b-y) 

(A-y)^/7r^ 


/W-y)(b-yj 
(x<- y)"*- /T-y  •' 


(xy-  A--M’)'/(Q.-y)(b-y)  ^ 
(A-y)*'*-  (/M-y)'''  Tr 


CLt-  b  ~2  y 

vTi  -  y )  ( b  -  yK  y  -  AH  y-x^o 


yy 


(1^)-  - 


(234) 


Conical  Coordinates 


(235) 


v^- 

X  = 

8"  ^  ^ 
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X'-- 


(236) 
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(237) 

(238) 

(239) 

Continuity  Equation. 

^  yj  + 

i  9  V 

A  M^- V'*  fV/4}-^ 

i  (/(b*"- V‘)(c’-*y‘)  ^  f - - -  ./  \ 

> "  <iy  Um'--  fVy)-  0 

(240) 

Energy  Equation .  The  left  side  of  equation  (5) 
becomes 
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-/w‘  I/'  l/»^l 4.  ^ 3.  r,  ^ 

Z (  Va  < IjM  +  l/y  )J  J  1  ^'•-y-  a/<  L  tn  + 

-  L/ >  V‘  1/0/  J.  ](i>'-»'>(c'-y»  3.  f , 

’’’ II /i  i  M'"  >''■  3>>Ln* 


(241) 


Also, 


v,f? ")]  -  (/a 

Vy/T  *''')]  +  ^  [if/M’-'V*  (v'AfT''*’+ 


^Kf.  =  ^  [VaFa  +  V'a.Fa,  +  /yFy]  , 


(242) 


(243) 
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The  components  of  the  stress  tensor  are 
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+ 


i/a'-y')Cc'-y»-9  f, - - 


(2ne) 


tr 


MM 


(>2-  "ir)[j-lA(A'-\;;)'^ 


ACM‘-y*)  3/1  U/<-y‘' 


/(b^-y*)^g^-y>-)  ^  / , _ , 


2m  _ 

A  UM^-bMCc'-AiMCA^'-y*)  Km  “ 


zy  /c6*-yo(c»--y^) 


(249) 


and 
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iq-  it)  [  A"  <iA  ■♦“ 


/  (m^-  b*Kc^-At^)  ^  ^ 

A(x'-v‘) 


A^M^-y“o  ay(C^v;)j  + 

r  2  ^  /  V  U  -  1/  4. 

vl  A  ^yW  x»--y*-  V>;  j  A  •♦■ 

2^  b^  4  c^-  ~  z  y^  ^  ■] 

A  V(b'-y')(c^-y')(x‘'y')  '^yj  . 


(250) 


Equation  of  Motion .  The  terms  of  Interest  In  the 
equation  of  motion  are 


^Va 

at 


t 


^Va 


4- 


h.  n^Zalli^lzJLD  V' V'/ 

A  A<-"-y"  ay  “a  a  ^  (251) 
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'ii,K  "  AV  ><'•- 3a  A  (a  IV'”  )  -^ 

^  Rb^y^SE?^^  />^«>/\  A/S^i27c5i?r . A 

b‘J  ( C  '-><>0  (A  >>>3  - 

A  (;u^ tr  ^  A  ^  “ 


2ji^E527cn?  ^  .y 


(255) 


■4=^  A  i  fib'--v)(c'^-y)  3  p  I  9  /  A  ji\ 

Ui.K  =  dS;  1- A  aAlx^  '^j  4- 


I II^^)  l.  ^ 

13# 

3i  U(b'-»')(c'--y'  'fr^")  -1- 
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A  report  of  this  nature  has  no  conclusions  in  the 
usual  sense  of  the  word.  The  methods  used  in  arriving 
at  the  final  results  were  not  new.  It  was  known  at  the 
outset  that  they  would  produce  valid  rf^siPts  if  they  wre 
applied  correctly.  This  project  was  undertaken  because 
the  equations  of  fluid  mechanics  are  not  readily  avail¬ 
able  except  in  terms  of  three  or  four  oi‘  the  more  common"' 
used  coordinate  systems.  The  tabulated  equations  in  Sec¬ 
tion  III  are,  therefore,  real"','/  the  "conclusions"  of  this 
thesis . 
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Appendix  A 

Description  of  Coordinate  Systems 
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The  coordinate  transformations  presented  here  were 
taken,  with  minor  changes,  from  the  work  of  Chester  n.  Pare 
(Ref  3:  93-107). 

( 1)  Cartesian  coordinates 

X  »  X  -00  <  X  «  00 

y  e  y  -00  i  y  i  OO 

Z  -  Z  -00  <  z  t  CO 

Sn  *  Svv*  (jJJ*  - 

and  the  surfaces  x.  y,  or  z  =  constant  are  planes. 

(2)  Cylindrical  polar  coord Inates 

X  -  r  cos  0 
y  =  r  sin  © 
z  s  z 

with  the  variable  ranees 

-OO  s  z  s  oo 
0  ^  r  £  00 

0  f  0  i  fr. 

The  surfaces  are: 


z 

=  cons' 

rsj 

p'  r?nc3 

r 

=  const 

cy i ' no  era 

e 

=  const 

■^’rnes  '  lirour': 

t'li  “  '.ijj"  j  • 


Wo  have: 


GAEA:ii/62-4 
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(3)  Spherical  polar  coordinates 

X  =  r  aln  e  cos  qp 
y  =  r  Bind  slny 
z  =  r  cos  ^ 

0  s  r  <  oo  ,  r  r  const  spheres 

0  s  </>  <  2V,  (f  =  const ~  azimuthal  planes 

0  5  0  -  )r,  6  -  const  circular  cones 

9..--  i 

(4)  Parabolic  coordinates 

Two  sets  of  coordinate  surfaces  are  generated  by 
rotating  the  parabolas  of  Pig.  1  about  the  x  axis  which 
Is  then  renamed  the  z  axis.  The  third  set  of  coordinate 
surfaces  are  azimuthal  planes  through  the  newly  relabeled 
z  axis. 

The  transformations  are: 

X  -  /XJk  cos  ^ 


y  =  /Xm  sin  <p 


0  <  A  £  oo  j 

0  <  oo  j 

0^ 


A  -  const  /V-  paraboloids: 
x^'^-  y  +  2Az  =  a’' 

=  const  paraboloids: 

x'^+  y  -  2^z  = 

(p  -  const  azimuthal  planes 


9n  = 


A^^M 


gjj  =  A/h 
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(5)  Prolate  spheroidal  coordinates 
This  coordinate  system  Is  constructed  by  rotating  the 
curves  of  Pig.  2  about  the  x  axis  which  Is  then  renamed 
the  z  axis.  The  coordinate  surfaces  which  are  generated 
are  prolate  spheroids  and  hyperboloids  of  two  sheets.  The 
third  set  of  coordinate  surfaces  are  azimuthal  planes 
through  the  nev;ly  renamed  z  axis. 

If  we  let  A  =  cosh  ,  /<=  cos  /S  ,  we  have 
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A<  oo, 

1  j 

o  <  <  2,tr, 


X  =  a/ >'■-  1  1/ 1  -  m'"  cos  4> 
y  =  a/A'"-  1  fi  -  /Ay  sin  9 


z  r  eAm 


A  =  const prolate  spheroids: 


I."  - 

A** 


a*- 


>(  const hyperlolds  of  two  sheets: 


9  -  const  /\^  azimuthal  planes 


S« 


Sxv-- 


J 


3„=  a.’-(A'- 


(6)  Spheroidal  coordinates  (oblate  spheroids) 

This  coordinate  system  is  constructed  by  rotating  the 
curves  of  Pig.  2  about  the  y  axis  which  is  then  renamed 
the  z  axis.  The  generated  coordinate  surfaces  are  oblate 
spheroids  and  hyperboloids  of  one  sheet.  The  third  set 
of  coordinate  surfaces  are  azimuthal  planes  through  the 
new  z  axis. 

X  =aA^  cos  (p  or  r  =  a  coshP^cos/? 
y  =  a  A  sin  (p  z  =  a  slnhev^sin^ 
z  =  a 

A  =  const oblate  spheroids: 


A*--! 


=  CL* 
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(7)  Parabolic  cylinder  coordinates 
The  traces  of  the  coordinate  surfaces  of  this  coordi¬ 
nate  system  on  the  xy  plane  are  shown  in  Fii^,  1, 

A  *  ^ 

X  =  X 


z  -  z 

0^  X?  a>j  A  =  const  /v  parabolic  cylinders,  ZAX  -  A'’ 
OSM-COj  M-  const/w  parabolic  cylinders,  lj'- 
-0o£2seD,  2  "  const planes 


a+_m 

Sm= 


A  t  ^ 

> 


GAE/M"/62-n 


The  traces  of  the  coordinate  surfaces  of  this  coordi¬ 
nate  system  on  the  xy  plane  are  shown  in  Fig.  2. 

X  c  a  /a'-J or  x  ?r  a  sinh®<8in  ^ 

y^aAM  y»a  coshYcos  ytf 

z  r  z  z  =  z 

.  _X'  , 

A  =  const  /V  elliptic  cylinders,  yt.  *  A 


M 


consi /'v- hyperbolic  cylinders, 


»  dL*- 


7  -  constA>  planes 
■  with  the  variable  ranges 

-O0S€y'-+oOj  1  ^  X  - 

-Od  Sf  i  ^  flo  . 


We  have 


d.'- 


1 


(9)  ellipsoidal  coordinates 

^  =  feL'-- c,')  Uin( 


-to<Xf  Cj,<jn< 


1. 
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The  surfacoa  are: 


A  "  const  /V  elllcsoid: 


Ja2  . 

^  b'-X  t^-A  "ij  a.>^)>C,-a)£A^c''j 


M=  const /v»  hi'pcrboloid  of  one  sheet: 


b'-  2’-  . 


6^-/4  "**  ^  ^  -  /K  - 


y  -  const  hyperboloid  of  two  sheets: 

X  ^  'b'*  ? ' 

0^  +  cTh-  ’  -i  ,  /)■■£  Vs  a.’ 


Also, 


i  ^Ai-A)(y-A:) 

^  (a'-'A^Cb^-A^Cc^A)  5 

i  (A-/mU>^->m) 


(lO)  Confoca 1  naraboilc  coordinates 


GAE/r^S/62-4 

The  transformations  are: 

1  , 

>=  Z  (XfM  t  y  -  d.-^) 

4  '  ■ 

3  2  Cj>~AH  b-At)  (  b  •  v) 

^  ^  a  -  b  . 

The  surfaces  X,  M>  and  constant  arc  the  paraboloids 

U  a. 

A  A(A-dL)  M  A  -  b)  =  1 

ZJf  ,  ij"*-  5ft 

u  +  -iL.  < 

y  y(v-4,)  Jh^b)  ’ 

with 

“00  ^>^<G-<M<b<A«=Oo, 

Also, 

^  .  A  (A-^)a-y^ 
ij,,  -  4  (a.-AXh-.O  , 
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i  (m-^)  (M-if) 
^J.1  “  ^  (^~M)(b-M)  ) 


anri 


1  (>^-A^(y-i4) 

’  4  (0. - yX b - . 

( 1 1 )  Conlca  1  coorciJ  nates 

The  surfaces  A,  m,  and  constant  arc  the  spheres 
x"*"  +  +  z  =  A*"  and  the  two  cones 


X'-  2-^ 

+■  yvT^’'  ■'■  ""  0  ^ 

whore  o^y<bs>L^<:C  J  O  <  CO. 

The  conoSj  constant.  Intersect  the  planes,  z  =  constant. 

In  ellipses.  The  cones,  V-  constant,  intersect  those  z 
planes  in  hyperbolas,  but  intersect  the  planes,  x  «  cons¬ 
tant,  in  ellipses.  Hence,  wc  can  visualize  the  cones  as 
being  elliptica]  cones  centered  about  the  z  and  x  axes. 

The  transformations  are 

X  =■  ) 
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Appendix  3 

Derivation  of  Basic  Equations 
of  Fluid  Mechanics 
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The  derivations  shown  here  can  be  found,  In  slightly 
varying  form,  In  the  works  of  Sokolnlkoff  (Ref  4:  290-324) 

and  McConnell  (Ref  1:  270-282),  among  others.  These 

specific  forms  were,  however,  flrsu  shown  to  the  author 
by  Lt.  Ray  M.  Bowen  of  the  Mechanical  Engineering  Depart¬ 
ment  of  the  Institute  of  Technology. 

Equation  of 'Continuity 

The  equation  of  continuity  can  be  expressed  in  inte¬ 
gral  form  as  follows: 

=  0  .  (B.I) 


Application  of  Leibnitz's  Rule  yields 


(B.2) 


By  the  Divergence  Theorem 


(B.3) 


Then,  since  volume  is  arbitrary, 

ll  +  =  0, 


(B.4) 


but 
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(  dx‘(<^  f  .  (B.5) 


(Ref  1:  155).  The  continuity  equation  then  becomes 


1/  +  q  akdq 


(B.6) 


The  physical  componcn’.s  of  the  contravarlant  quantity  tT  ‘ 
are  represented  by  f^am)  If  we  denote  the  physical  com- 
poncnts  by  the  symbol  v(  ,  then  cr  « then  write 
the  continuity  equation  in  its  final  form: 


It  +  /g  ax‘  f  ^  • 


(B.7) 


Energy  Equation 

The  energy  equation  can  be  written  in  integral  form 
as  follows: 

u- ir‘ 'l  o'!-  =  4/)  J 

'j  ^  d't  f  of  5 .  (b.8) 

Application  of  Leibnitz's  Rule  yields 
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lai  dY  +  I,,,  [f{uf  i  t/;lf‘)i/hjds= 

dst  ds^  (B,9) 

and,  by  the  Divergence  Theorem, 

X  v.i;*)] +  ff  (ut  f  Lr.tr‘V‘],j}=c('K 

X«)  f  ^  (B.IO) 

Then,  since  the  volume  is  arbitrary, 

If  [f  (UtflTiU-')]  t  [j>(ut  = 

+  fir,  F‘~  .  (B.u) 


By  expanding  and  subtracting  out  the  continuity  equation 
we  obtain 

f  [lt(wf  iu-u')  +  V*(u+i(riV‘),*J  ■- 

*■  ftfi  F'  -  .  (B.12) 

But 
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U  - 


(B.13) 


Substitution  of  (B.13)  into  (B.12)  yields 

J  =  -  (pU-^),  ;  + 

f^iF‘  -  ,  (B,i4) 

v;here  v/e  have  also  made  use  of  the  relationship 

^  =  -  p  .  (B.15) 


After  simplifying  and  rearranging,  the  energy  equation 
becomes 


I?  +  +  f  ,  (b.16) 


or,  in  terms  of  physical  components 
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+  I  V'iK)+  75^  *iKK)]  -  it  *■ 


^  ix‘  Ki  )r  ‘■^ )  ^  f  Fi  - 


V*M 


!(3J7) 


:£quallon  of  Motion 

The  equation  of  motion  can  be  written  in  integral 
form  as  follows: 


(B.IS) 


where 


4.  . 


(B.15) 


By  Leibnitz's  Rule 


Jyrt,  It  dt-  = 


(B,.19) 


andj  by  the  Divergence  Theorem, 
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/y,„  f  ft  (  f  1/  ‘  )  +  (  If  d  t-  '  %  (B.  so) 

Then,  since  the  volume  is  arbitrary, 

at  (fv‘)+  (^v‘vh,j=  »•  jF\  (B,2i) 

Ho’w,  after  expanding  and  subtracting  out  the  continuity 
equation,  we  have 

or 

f  -f  -  f  ^‘  =  'h'i  j  >  (3.23) 

where 

i.  -  at  ^  ,  j  e  =  jt  .  (B.2H) 

Equation  (B.23)  can  be  expressed  in  covariant  form  as 

(B.25) 
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Appendix  C 

The  Three  Basic  Equations 


94 


GAE/ME/62-4 
Continuity  Equation 


Energy  Equation 


f  ViV;)+ *('(/?+ i  =  H  + 


Equation  of  Motion 


(12) 
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